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O . ABSTRACT 

Context. Radiation driven winds are the likely origin of AGN outflows, and are believed to be a fundamental component of the inner 
structure of AGNs. Several hydrodynamical models have been developed, showing that these winds can be effectively launched from 
AGN accretion discs. 

Aims. Here we want to study the acceleration phase of line-driven winds in AGNs, in order to examine the physical conditions for the 
existence of such winds for a wide variety of initial conditions. 
Methods. We built a simple and fast non-hydrodynamic model QWIND, where we assume that a wind is launched from the accretion 
disc at supersonic velocities of the order of a few 10 2 km/s and we concentrate on the subsequent supersonic phase, when the wind is 
. accelerated to final velocities up to 10 4 km/s. 

Results. We show that, with a set of initial parameters in agreement with observations in AGNs, this model can produce a wind with 
terminal velocities of the order of 10 4 km/s. There are three zones in the wind, only the middle one of which can launch a wind: in 
the inner zone the wind is too ionized and so experiences only the Compton radiation force which is not effective in accelerating gas. 
This inner failed wind however plays an important role in shielding the next zone, lowering the ionization parameter there. In the 
middle zone the lower ionization of the gas leads to a much larger radiation force and the gas achieves escape velocity This middle 
zone is quite thin (about 100 gravitational radii). The outer, third, zone is shielded from the UV radiation by the central wind zone 
and so does not achieve a high enough acceleration to reach escape velocity. We also describe a simple analytic approximation of our 
model, based on neglecting the effects of gravity during the acceleration phase. This analytic approach is in agreement with the results 
of the numerical code, and is a powerful way to check whether a radiation driven wind can be accelerated with a given set of initial 
OO ' parameters. 

, Conclusions. Our analytical analysis and the fast QWIND model are in agreement with more complex hydrodynamical models, and 

allow an exploration of the dependence of the wind properties for a wide set of initial parameters: black hole mass, Eddington ratio, 
f*"^ ■ initial density profile, X-ray to UV ratio. 

Key words, quasars - winds - models - X-rays - broad emission lines. 

1 . Introduction applied to accretion discs in compact binaries ( Proga, Drew & 

Stone 1998) and to AGN discs (Murray et al. 1995, Proga et al. 
Outflowing winds are now believed to be common, and quite 200a hereafter p 0a Proga 2003) xhese authors use the power . 

possibly ubiquitous, in the inner parts of Active Galactic Nuclei M hydrodynamical code ZEUS2D (Stone & Norman 1992) to 

H , (AGNs) and quasars. The most striking evidence comes from solye ±e wind ations> assuming a Shakura-Sunyaev (1972, 

the 10% of Broad Absorption Line (BAL) quasars which show hereafter ss) a . disc> and plausible initial conditions. The main 

broad absorption lines blueshifts spanning 10-20 thousands resuhs of thefje works m (1) ±e demonstration that a wind 

km s- 1 . Less spectacular, but more common, evidence for out- C£m fee launched and acce i e rated up to velocities of the order 

flows comes from the 50% of AGN with Narrow Absorption of 1q4 km ^ and (2) the determination of the wind geometry 

Lines (NALs) both UV and X-ray (Reynolds 1997, George et al. and h kal ^ for ±e iyen starti conditions . These mod . 

1998, Crenshaw et al. 1999, Krongold et al. 2003, Vestergaard dg haye been tested wkh seyeral different choicefj of ^ 

2003, Piconcelh et al. 2005, Ganguly & Brotherton 2007). cond itions (Proga & Kallman 2004), showing that a wind can 

Several attempts have been made to explain the origin of ^ for a wide of black hole masfj and accretion rates 
such winds. Two main scenarios have been suggested: magnet- 
ically driven winds (Blandford & Payne 1982, Konigl & Payne Further recent improvements in this field are the modeling 
1994, Konigl & Kartje 1994, Everett & Murray 2007), or ra- of Schurch & Done (2007), where the interaction between the 
diation driven winds. One of the most promising explanations X " ra y radiation and the outflowing wind is throughly analyzed, 
is through radiation line-driven winds arising from the accre- and the anal y sls of possible effects of radiation winds at larger 
tion disc. The physics of radiation line-driven winds based on distances from the accretion disc, such as on a parsec-scale X- 
the work of Sobolev (1960) was developed by Castor, Abbott & ra y heated toms (Dorodmtsyn et al. 2008), and on large-scale 
Klein (1975, hereafter CAK), and Abbott (1982, 1986)forwinds AGN outflows (Kurosawa & Proga 2009). 
from hot stars. CAK showed that resonance line absorption in an One key aspect of the physics of radiation driven winds, 
accelerating flow could be hundreds of times more effective than clearly emerging from the models mentioned above, is that re- 
pure electron scattering. More recently the same theory has been gardless of the details of the initial launching phase, it is im- 
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possible to avoid the strong, probably dominant, effect of radi- 
ation pressure in the subsequent acceleration phase, where the 
wind gains more than 99% of its kinetic energy. This is easily 
estimated from the comparison between the amount of momen- 
tum absorbed by the gas and its final momentum. For example, 
Hamann (1998) estimates that at least ~ 25% of the UV radiation 
emitted by the Broad Absorption Line Quasar PG 1254+047 is 
absorbed by a gas with column density of the order of 10 23 cm" 2 . 
It is sufficient that the luminosity of this source is of the order of 
10% of the Eddington luminosity to conclude that the momen- 
tum in the outflowing wind is of the order of that absorbed in the 
UV wavelength range. 

The distinction between the launching and acceleration 
phases of an accretion disc wind, which we made above, 
allows the modelling of the problem to be separated into 
these two parts. The launching phase has been investigated 
in several numerical simulations, both in hydrodynamical 
(Ohsuga et al. 2005) and magneto-hydrodynamical (Hawley 
& Krolik 2006) regimes. Recently, the creation of out- 
flows from accretion discs has been investigated through 
MHD simulations including the effects of radiation pressures 
(Ohsuga et al. 2009). 

The main aim of the work presented here is to explore a wide 
set of initial conditions for the acceleration phase of a wind in 
AGN using a deliberately simplified approach in the hope that 
this can produce an intuitive understanding of quasar winds and 
so guide future detailed simulations. We will make use both of 
analytic approximations and of a simplified numerical code for 
quasar winds (which we named QWIND). QWIND treats in de- 
tail the radiation force mechanism, but not the internal gas pres- 
sure in the wind. As the wind velocity in the acceleration phase 
is always many times the thermal velocity of the wind gas, this 
is a reasonable approach. 

Our study is motivated by the observations of fast outflows 
in Broad absorption Line (BAL) quasars (Weymann 1997), by 
photometric evidence for a highly flattened structure in BAL 
winds (Ogle et al. 1999) and by results suggesting an axially 
symmetric, but not spherical, spatial distribution of the Broad 
Emission Line (BEL) gas (Wills & Browne 1986, Brotherton 
1996, Maiolino et al. 2001, Rokaki et al. 2003). Elvis (2000) 
has proposed a specific quasar unification model based on this 
kind of structure. This model predicts that a geometrically thin, 
funnel-shaped outflow is present in all AGNs and quasars, which 
allows the observational properties of the different classes of 
sources to be explained through orientation effects. This model 
provided our initial motivation, as we suspected that a radiation 
driven wind might naturally create a thin, funnel-shaped struc- 
ture. 

We therefore developed the QWIND model in order to easily 
and quickly explore the acceleration of a radiation-driven wind 
for the huge range of physically possible initial conditions, in 
terms of black hole mass, accretion rate, relative strength of the 
X-ray radiation, initial density and temperature of the outflowing 
gas. Our approach does not add more physical insights for the 
single wind solutions, than the already available hydrodynamical 
models mentioned above, but allows a much wider exploration 
of the initial parameters space. 

The purpose of this paper is to present the model code, 
QWIND, and to demonstrate that QWIND produces results in 
agreement with both an analytic treatment and with the more 
complex and detailed approach of POO. The structure of this pa- 
per is the following: in Section 2 we briefly review previous 
results on the topic of radiation-driven winds. In Section 3 we 
present the QWIND code and we show examples of its possible 



applications. In Section 4 we present an analytical treatment of 
the wind equations, which provides interesting constraints on the 
existence of such winds, as a function of the initial parameters. 
In Section 5 we briefly compare our results with those obtained 
by more complex hydrodynamical codes. Finally, in Section 6 
we present our conclusions and outline future work. 

2. Summary of previous results 

The radiation force on a moderately ionized gas due to incident 
ultraviolet radiation is primarily via line absorption, rather than 
continuum electron scattering (CAK). In a constant velocity gas 
the 1000 times larger cross-sections of resonant line transitions 
over continuum Thomson scattering has little effect because the 
narrow wavelength ranges spanned by each transition contain 
relatively little UV continuum. However, a real wind, driven by 
a central radiation source will be accelerated, and so the wave- 
length for each absorption line will be Doppler shifted by the 
relative velocity of the gas with respect to the source. If the 
radial velocity gradient shifts the absorption wavelength suffi- 
ciently for each gas element, then fresh UV continuum is ab- 
sorbed, leading to continued acceleration. As the process contin- 
ues the fraction of the UV continuum that is absorbed is greatly 
increased. The absorption cross-section then becomes a local 
function of the physical conditions of the wind (the "Sobolev 
approximation", Sobolev 1960). CAK introduced a numerical 
"force multiplier", M(t), which represents the enhancement of 
the radiation force due to line absorption with respect to pure 
electron scattering. M(f) depends on only one local quantity, the 
"effective optical depth", f, defined as: 

dv r i 

t = o- T nv th \ — | (1) 
ar 

where crj is the Thomson cross section, n is the gas density 
v t h is the thermal velocity in the gas, and v r is the radial compo- 
nent of the velocity. 

Taking into account the thousands of, mostly weak, UV ab- 
sorption lines for each ion, CAK estimated the following ana- 
lytic power-law approximation for M(t): 

M(t) = Kxr a (2) 

where a is a positive number (with a typical value a = 0.6) 
and K is of the order of 0.1-0.4. 

However, as a gas becomes more ionized, either collisionally 
due to an increased temperature (e.g. through photoionization 
due to the intense UV source), the number of transitions avail- 
able for absorption is reduced quite dramatically. In launching 
AGN winds this is a well known problem, and ad hoc schemes 
to shield the gas have been proposed (e.g. "hitchhiking gas", 
Murray & Chiang 1997) 

A useful, second-order approximation, which takes into ac- 
count the ionization state of the gas and solves the divergence of 
the previous equation for small f, is the following (Abbott 1986): 

M(0^xr-x[ (1+T ^r ? " 1 ] (3) 

T MAX 

with tmax = t x tjmax, where t\max is a coefficient that is 
strongly dependent on the ionization factor, This relation was 
studied by Stevens & Kallman (1986) who produced an analytic 
approximation that is correct within 10% over the whole range 
of interest for We plot this relation in Fig. Q] together with 
the dependence on f of the constant K. In a wind, £ can be as 
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Fig. 1. Analytic approximation fo/the parameters K (solid line) 
and r/MAX (dashed line), as a function of £ (Stevens & Kallman, 
1986) 
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Fig. 2. Force multiplier M(t) vs optical depth t for a neutral gas in 
the Castor, Abbott & Klein (1975) formulation (straight dashed 
line), and the Abbott (1982) and Stevens & Kallman (1986) for- 
mulation that allows for ionization state and is more accurate for 
low effective optical depth (t < 10~ 5 ). 



high as 10 3 - 10 4 in the inner surface, but then rapidly drops 
to small values (£ < 1) because of the lowered X-ray flux. As a 
consequence, in the range of interest for a wind both t/max and K 
are constant to within a factor ~ 2. In Fig. [2] we plot Eq. 2 and 3 
for a neutral gas. We note that the correction factor is significant 
only if t < 10" 5 . 

Using these results, POO developed a code which solves the 
hydrodynamical equations for a wind arising from an accretion 



disc, pushed by a central radiation force described by the above 
equations. In this work, a SS disc is assumed, and the ionization 
state is determined by the incident X-ray radiation from the cen- 
tral source. The main conclusion of POO (in agreement with the 
previous work of Murray et al. 1995) is that a wind can easily 
arise from the accretion disc of AGN, under some initial con- 
ditions (SS disc, with a choice of free parameters such as the 
Eddington ratio and the X/UV ratio, which will be discussed in 
the next Sections). The dependence of the wind properties on 
the black hole mass and accretion efficiency is particularly in- 
teresting for they can now be tested against observations, as we 
discuss below. The black hole mass, Mbh and the accretion effi- 
ciency relative to Eddington, €edd, are the parameters determin- 
ing the luminosity of AGN. 



3. Radiation-driven acceleration of a wind 

Several physical effects can be considered concerning the equa- 
tions of the radiation force, to understand whether this exter- 
nal force is capable of accelerating the wind up to velocities of 
~ 10 4 km s or higher, as observed in BAL quasars. 

The results summarized above (Section 2) show that a ra- 
diation driven wind can arise from accretion discs in AGNs. 
However, the breadth of initial conditions that produce such a 
wind are at present unknown. First, while the SS disc is a plau- 
sible solution, many other disc structures are possible. Secondly 
and most importantly, even in the SS disc paradigm several ini- 
tial parameters are unknown: the black hole mass, the Eddington 
ratio, and the viscosity parameter a. Finally, since the physical 
reason for the initial launching of a wind is not clear, the radius at 
which the wind is launched is unknown. However, recent stud- 
ies of warm absorber winds in local AGNs suggest that winds 
on scales as small as a few thousand gravitational radii from the 
central black holes are possible (Krongold et al. 2006). 

An important aspect of the solutions found in the previous 
works of Murray et al. (1995) and Proga et al. (2000) is that 
the sonic and critical points of the wind are reached at a small 
height H above the disc with respect to the distance from the 
central source R (H/R < 0.1). Also, the velocitiies at the sonic 
and critical points are of the order of a few 10 2 km s _1 , which 
is only a few percent of the final BAL velocities, of the order of 
lf/kms -1 . 

As briefly discussed in the introduction, radiation pressure is 
expected to be relevant (and probably dominant) in the wind ac- 
celeration, regardless of the details of the launching mechanism. 
However, following the theory summarized in Section 2, in or- 
der for the absorption line mechanism to be effective, a balance 
between the intensity of the X-ray radiation (which determines 
the ionization parameter) and the intensity of the UV radiation 
(which provides the momentum needed by the gas to accelerate 
to high velocities) is required (Murray & Cheng 1995). 

In this paper we discuss the general conditions for the 
acceleration of a line-driven wind after the launching phase. 
In particular we investigate whether the "thin wind" geometry 
proposed by Elvis (2000), which is successful in explaining a 
wide set of observational properties, can arise naturally. The 
main unexpected features of the Elvis (2000) geometry are that 
(1) the wind is thin (AR/R<1) and that (2) the wind is initially 
quasi-vertical, making a hollow cylinder of height h/r ~ 1, 
before becoming a radial, biconical flow. 
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3.1. The QWIND code 

We developed the QWIND code with the following assumptions: 

1. The X-ray source is point-like, isotropic and located at the 
center of the disc. 

2. The UV/optical source is the accretion disc, emitting ac- 
cording to the SS model. In our model the wind inner 
radius can be as small as 100 R s , so the dimensions of 
the emitting region of the disc are non-negligible, and 
the disc cannot be considered a point source. Therefore, 
the dependence of the disc emissivity with radius, and 
the actual disc extension are taken into account. We also 
take into account the anisotropy of the disc emission, 
which has important consequences for the wind proper- 
ties, the radiation being zero on the disc plane, and in- 
creasing with the cosine of the inclination angle of the 
disc with respect to the wind. Limb darkening (Fukue 
& Akizuki 2007), which would accentuate this effect, has 
been neglected. 

3. The gas arises from the disc at all radii between an inner 
radius R/n and an outer radius Rout, at a height H above 
the disc, with a given density «, temperature T, and initial 
vertical velocity Vq. The radial velocity of this gas is zero, 
while the angular velocity in the disc plane is Keplerian. 

4. The wind is subject only to the central gravitational force 
Fq due to the black hole of mass Mbh and radiation force 
F r . Internal gas pressure is neglected, as is the gravity of the 
disc. 

5. The radiation force in each element of the gas is computed 
using the CAK and Abbott (1986) equations described in 
Section 2 (Eqs. 2 and 3). The ionization parameter, is de- 
termined assuming a cross section <x x = 100<x r if £ < 10 5 , 
and ctx - cr T if £ > 10 5 . This simple approximation re- 
flects the fact that when the ionization factor is too high, 
there are no bound electrons for photoelectric absorption to 
be effective. Below this critical value, photoelectric absorp- 
tion on the electrons in the inner shells of metals is much 
more effective than Thomson scattering in removing X-ray 
photons. Using this approximation, if £ > 10 5 , X-rays pene- 
trate deep into the wind, keeping £ high, until the Thomson 
optical depth is higher than 1 (Nh > 1.5 x 10 24 crrT 2 ) or the 
R 2 factor makes £ decrease below 10 5 . Until this point, the 
radiation force is negligible because M(t) is small. Deeper 
inside the wind, where Tx > 1, the UV radiation is also fully 
absorbed, and again <x x = crj-, and therefore the radiation 
force will be negligible in this case too. If instead £ < 10 5 
the X-rays are absorbed in a thin layer Q(JVh ~ 10 22 cm" 2 ) 
and the ionization parameter drops as £ oc e -iooa- T N H ^ la p[^\y 
reaching the values at which the force multiplier becomes 
significantly higher than 1 . 

QWIND follows a gas element of radial thickness AR (A 
R<<R) from the accretion disc to radii greater than 10 3 Rs , start- 
ing from a set of initial radii This approach, consisting 

1 The problem can be complicated by the effects of secondary ioniza- 
tions, which make f decrease more slowly than exponentially (Maloney 
et al. 1996). However, this effect is not important when f is high, and 
the external electrons (the targets for secondary ionizations) are already 
free. 

2 Note that throughout the paper Rj and n - ( refer to the initial radius 
and density, respectively, of the i-th stream line, while R 1N and n 1N are 
the initial radius and density of the first stream line, i.e. Rj N =Ri(i=l), 
w ;jv =n,(i=l) 



of following the trajectories of a fluid element, is similar to 
that adopted in several previous works analyzing disc out- 
flows (Icke et al. 1980, Tajima & Fukue 1996, Watarai & 
Fukue 1999). 

The acceleration of a gas element belonging to the i- 
th stream line, at coordinates R, <f>, z, is simply given by: 
a = -GM BH R~ 2 + F RAD /m(R) where m{R) is the mass of the 
gas element, calculated from the mass conservation equation 
(therefore depending on the radius and velocity of the gas el- 
ement, and on the initial density profile). The radiation force 
Frad is estimated integrating the contribution of each disc 
element: 

Frad = f f(cf>)e- T -^—o- T (M(t) + 1) (4) 
J disc c 

where f(<p) is the geometrical projection factor of each disc 
element as seen from the gas element; / is the luminosity of 
the disc element; M(t) is the force multiplier, estimated as 
described above, and r is the optical depth, estimated adding 
the contributions of the stream lines between the gas element 
and the disc: T=E,cr 1 xAR,xn(/?,)/f(0), where n(/?,-) is the den- 
sity of the gas element in the i-th stream. The luminosity / 
of each disc element depends on the radial distance from the 
center of the disc, as estimated in a SS disc. 

Our choice to integrate over the disc structure is moti- 
vated by the non-negligible dimensions of the disc with re- 
spect to the inner wind streams: in our computation the 
disc emission is calculated up to an outer radius of 400 R s 
(even if its contribution to the total luminosity is negligible 
at R>100 R s ), while the inner wind radius is 100 R s (see 
next Section for further details on the computation domain). 

The computation starts from the inner stream line (free 
from absorption), and then proceeds to the outer ones. At 
each step, the calculations are done in the following order: 

- We estimate the absorption due to the inner stream lines. 
For the UV component, this is done simply adding the col- 
umn densities of the inner stream lines, and multiplying by 
the Thomson cross section. For the X-ray component, the 
computation is more complex: for each stream line crossed 
by the light ray from the disc element, we estimate the ion- 
ization factor and then, depending on its value, the X-ray 
cross section, following the prescription described above. 

- Using the absorption correction estimated in the previous 
point, and integrating over each disc element, we compute 
the flux on our wind element. 

- We estimate the density of the gas element from the val- 
ues of density, mass and velocity of the previous computation 
step, and requiring mass conservation. 

- We estimate the force multiplier M(t), following the pre- 
scription described in the 5th point of the code properties. 

- We obtain the acceleration on the gas element, and we com- 
pute the new spatial coordinates and velocity. 

As a consistency check of the numerical stability of 
QWIND, we performed several runs with no radiation force 
contribution, and we checked that the evolution of the stream 
lines are those expected in a Keplerian problem. In particu- 
lar, we checked that the conservation of angular momentum 
holds, and that the results are stable against changes in the 
number and dimensions of the gas elements used in the com- 
putation. 

Neglecting the internal gas pressure makes the stream lines 
of gas in QWIND intersect unphysically. When two stream lines 
intersect, this means that the gas in the internal line is pushed 
outwards more than the gas on the external line. In reality this 
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obviously cannot happen, and is prevented by the horizontal 
component of the gas pressure. However, the two stream lines 
do delineate a limiting cone that a fully modelled wind must 
keep within. The gas pressure force between the two "lines" will 
make the gas move along a line whose inclination angle with 
respect to the disc axis is intermediate between the two lines 
computed by QWIND. This approximation therefore makes the 
results of QWIND usable only in this bounding sense, but should 
not change the general findings that we are interested in, i.e. 
whether the wind can be accelerated by radiation pressure, the 
approximate angles to the disc, and covering factor of the cen- 
tral source, since these properties strongly affect the observables. 
Another important effect neglected in our treatment is the change 
in momentum at the bending of the wind, which is expected to 
produce shocks and, possibly, soft X-ray emission, although this 
is unlikely to dominate the X-ray luminosity of the AGN. The 
case for the outflowing gas as a source of the soft X-ray excess 
observed in many quasars has been discussed by Pounds et al. 
(2003). 

The number of stream lines used in the runs discussed in 
this paper is kept to just 20. We also made several runs with a 
larger number of lines (50), to check for discretization effects 
and found that the properties of the solutions do not change sig- 
nificantly. 

The initial density of the wind at a radius R is one of the main 
free parameters of our model. In general, we parametrize the 
density as «,(/?,) = h/n x(/?,7/?/Ar)~' 3 . Two particularly interesting 
cases are (3 — (constant density at all radii) and [} = 15/8, the 
density scaling in the SS disc. In this paper we adopt f3 — in 
our examples. Different cases (and, in particular, the SS disc-like 
profile, will be investigated in a forthcoming paper. 

The initial vertical velocity vo is a free parameter of the code 
and is of the order of a few 100 km s -1 . Note that this velocity 
is much lower than the escape velocity, therefore most of the 
kinetic energy needed to make a wind must be provided by the 
external radiation force. 

With these assumptions, the total mass outflowing from the 
disc is easily derived, integrating the contribution from each disc 
ring. The fraction of this mass that falls back onto the disc ver- 
sus the fraction which escapes through a wind, depends on the 
subsequent acceleration phase, and is studied by our model. 

The initial parameters Vo and «,(!?/) are in reality a result of 
the previous launching phase. For example, in a purely radiation- 
driven wind, the mass loss rate is uniquely determined by the 
regularity and stability conditions at the 'critical point' (CAK, 
Lamers & Cassinelli 1999, Murray et al. 1995). Since we are not 
studying this launching phase here, a fundamental requirement 
is that our final results do not depend critically on the exact val- 
ues of these parameters. Since the density profile does affect the 
acceleration phase of the wind (because it determines, together 
with the X-ray flux, the ionization parameter), it is particularly 
important that the final results of the acceleration phase are in- 
dependent on the exact value of the initial velocity, provided that 
it is small compared with the final velocity, and large enough to 
be supersonic (this, for temperatures of a few 10 6 K implies Vo 
of the order of a few 10 2 km s _1 ). As we discuss in detail below, 
we have carefully checked that this is indeed the case for our 
wind solutions, and that, for particular values of the initial ve- 
locity (inside the range mentioned above), our solutions satisfy 
the critical point conditions for a purely radiation-driven wind. 

We next discuss the solutions obtained with our code 
QWIND for different choices of initial parameters. First, we con- 
centrate on a particular 'baseline' solution, in order to under- 
stand the physics of the wind. Then, in Section 3.3 we discuss 



Table 1. Baseline QWIND Model Parameters 



Parameter 


Value 


RlN 


100 R s 


Rout 


800 R s 


vo 


10 7 cm s -1 


£edd 


0.5 


M BH 


2x10 s Mo 


n, N 


2 x 10 s cm" 3 


P(n l (R) = n IN x(R,/R IN yl 3 ) 





T 


2xl0 6 K 


fx 


0.15 



the stability of our solutions. In Section 3.4 we discuss the de- 
pendence of the disc properties on the initial physical conditions. 
In Section 3.5 we briefly discuss a small survey of parameters in 
order to show the potential of our method in studying the wind 
properties in a variety of physical conditions. 

A systematic analysis of all the variables, including different 
density profiles, and the effect of toroidal magnetic fields, will 
be the subject of a forthcoming paper. 

3.2. Results: The baseline model, a case study 

We show here the results obtained from running QWIND with 
the initial parameters given in Table 1 . 

The inner and outer radii were chosen in order to investigate 
the region where a radiation-driven wind is expected, starting 
from an inner radius close to the outer UV-emitting region of the 
accretion disc, and studying the solutions up to an outer radius 
large enough to contain the inner broad line region. 

The initial velocity was arbitrarily chosen. We remind that 
the important points related to this parameter are: a) its value 
must be small compared with the escape velocity, and b) the ex- 
act initial value should not significantly affect the final results. 

A further fundamental parameter in the model is the ra- 
tio between ionizing radiation and bolometric emission, fx. 
Considering that the whole 0. 1-100 keV spectrum contributes to 
the gas ionization, and that the wavelength range of the intrinsic 
disc emission is from optical to X-rays, a typical value for this 
ratio is f x -15% (e.g. Elvis et al. 1994, Risaliti & Elvis 2004). 
Quite different values are however possible, due to the large dis- 
persion of the X-ray to optical ratio among quasars, and to its 
dependence on optical luminosity (e.g. Steffen et al. 2006). We 
will discuss the effects of changing this parameter in the next 
Sections. 

The main properties of the solutions are shown in Fig. [3] In 
the upper panel we show the wind lines in a plane perpendicu- 
lar to the accretion disc. In the lower panel we show the ratio 
between the final velocity and the escape velocity for each wind 
line. By definition, we have a wind when the gas velocity be- 
comes higher than the escape velocity, i.e. when the ratio in Fig. 
|3j) grows to values greater than 1. In Fig.|4]we show the velocity 
profile and force multiplier M(t) for the fifth stream line, i.e. the 
one reaching the highest final velocity. 

Several interesting properties of the wind can be drawn from 
Fig.|3]and|4] 

- Initial distance: Fig. |3j) shows that the baseline model 
gas reaches escape velocity only if it is launched between 
~ 1 50 Rs and ~250 Rs . At smaller radii, the gas is too ionized 
for the radiation force to be effective. At greater radii, the UV 
radiation is not enough to push the gas up to the escape ve- 
locity. The exact values for the allowed radius range depends 
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M 7 =10 n 7 =20 e E =0.5- 



f x = 15% 




Fig. 3. Results of the simulations assuming a constant initial 
density profile. In this example, Rjn=100 Rs, Rout=800 Rs, 
n /AI =2xl0 8 cm" 3 , v =10 2 km s -1 , /x=15%, e ED D=0.5, M 7 =20. 
The upper panel shows the wind geometry. In the lower panel we 
plot the ratio between the final velocity of each stream line and 
its escape velocity, as a function of the initial radius. We also 
show the final velocity versus radius for a wind with the same 
initial parameters but a larger inner radius (Ri N = 300 Rc) 



on the inner radius (in our case Rjn - 100 Rs ). If, for exam- 
ple, Rin = 300 Rs (and all the other initial parameters are the 
same), it would be impossible to have a wind at R = 300 Rs , 
because the gas is too ionized, due to the absence of shielding 
from the central X-ray emission. A wind at larger radii could 
however be possible, as we show in Fig. 3. QWIND thus pro- 
vides a simple explanation for the existence of a wind arising 
only from a small range of distances from the center. 

- Inner failed wind: the inner stream lines form a "failed wind" 
which shields the outer stream lines from the central X-ray 
radiation. This shielding is fundamental to decrease the ion- 
ization parameter in these stream lines, so allowing an effec- 
tive radiative acceleration. This inner component resembles 
the "hitchhiking gas" of Murray & Cheng (1995). 

- Radiation force: the upper panel of Fig. [4] shows the profile 
of the force multiplier M(t) for the fifth stream line in Fig. 5. 
Note that the maximum value is never higher than ~ 30. The 
relatively low values of M(t) (well below the values > 100 
estimated by CAK in hot stars) imply that only fairly large 
values of M/M^ (>0.03) can produce an escaping wind. 



> 




500 



1000 



5000 



Fig. 4. Force multiplier (upper plot) and velocity (lower plots) 
versus radius for the fifth stream line in our case study. The 
dashed line in the lower plot represents the escape velocity as 
a function of the distance from the central black hole. See text 
for details. 



A more technical consequence of the estimated values of 
M(f) is that the correction factor in Eq. 3 is never impor- 
tant, and Eq. 2 is a good approximation of the force multi- 
plier (Fig. 2, i.e. log t > -5). This makes possible an analytic 
analysis of the wind launching problem, which we discuss in 
Section 4. 

Velocity: Fig. [4] (lower panel) shows the velocity profile for 
the fifth stream line of our case study. The acceleration is 
fast, and the escape velocity is reached at ~ 500 Rs . Then the 
velocity continues to slightly increase. In this phase the ra- 
diation force is still effective in supporting the wind against 
gravity. Indeed, the "effective" Eddington ratio, which de- 
pends on the actual one, £edd, the force multiplier, and the 
inclination angle 6 of the disc as seen from the wind line, is 
e E DD x M(t) x cos 6 > 1 up to large radii (R~ 5000 R s ). 
Geometry: Fig. [3^ shows that the gas rises vertically for 
~ 50 Rs , and then bends toward a radial direction. Whether 
the gas falls down or maintains this direction depends on the 
effectiveness of the radiation force (see below). The angle 
of the wind above the disc is approximately 20 deg, giving 
a substantial covering factor of 35%. We note that the ac- 
tual value of the covering angle is somewhat uncertain be- 
cause the lines are treated as independent while, in reality, 
when two lines intersect, the gas pressure from the inner line 
pushes the gas on the outer line outwards, altering the final 
covering factor. 
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3.3. Check of initial conditions 

In order to test the reliability of our initial conditions, in partic- 
ular regarding the initial values of the vertical velocity and the 
height above the disc, we performed two checks: 

1 . We varied the initial height from H — 5 Rs to H — 15 Rs , and 
the initial velocity from 100 to 300 km s _I . No significant 
change in the final properties of the wind was found. 

2. We followed the wind streamlines backward in time down 
to the critical point (estimating this location as in Murray et 
al. 1995) and required that our solutions satisfy the regularity 
and stability conditions at the critical point (CAK). We iter- 
atively varied the initial velocity until a consistent solution 
was found. A solution was found with a value of the initial 
velocity within our range of study (vo=100-300 km s _1 ). 

The solution provided by QWIND thus represents the super- 
sonic part of a purely radiation-driven wind (i.e. the case where, 
even in the subcritical phase, the radiation force is the only ex- 
ternal force present in addition to gravity) and is in all respects 
equivalent to that found by Murray et al. 1995. The existence of 
such solutions within the range of initial velocities we are ex- 
ploring shows that our approach is well justified, at least in the 
pure radiation-driven scenario. 

3.4. Dependence on the individual parameters 

In order to explore the dependence of our solutions on the initial 
parameters, we proceed in two steps: First, we analyze several 
cases close to the baseline model, changing the density, the X- 
ray to UV ratio, the Eddington ratio, and the balck hole mass. 
Then, in the next Subsection we show the results of a survey 
of the £edd - Mj space for a few choices of the other initial 
parameters to show the potential of our method in exploring the 
initial conditions parameter space. 

First, starting from the baseline model (Table I3.2l i. we sys- 
tematically changed €edd, fx, n%, Mbh, one at a time: In Fig.|5][8] 
we show plots analogous to that in Fig.[3j5, i.e. the ratio between 
the final wind velocity and the escape velocity, v/vesc, against 
the initial wind radiusQFrom these results several additional in- 
dicators of the physical processes dominating the wind can be 
obtained, based on the dependence of v/vesc on the: 

- Eddington ratio £edd (Fig. 0): A too small value (panel A, 
£edd =0.1) implies that the radiation is inadequate to accel- 
erate the gas to escape velocity. With a higher value (panel 
B, £edd = 0.2, a stream line is able to slightly exceed the es- 
cape velocity, thus forming a wind. Finally, for much higher 
values (panel C, €edd = 0.9, to be compared with our case 
study with €edd = 0.5) a wind is effectively launched, but 
at larger radii than in the previous cases. This is due to the 
higher X-ray luminosity, which over-ionizes the gas up to 
larger distances from the center. 

It is interesting to note that the maximum velocity reached by 
the wind in the case eedd =0.9 is smaller than in the baseline 
case of Fig. [3] despite the higher Eddington ratio. This is due 
to the reduction of the UV radiation pressure on the stream 
lines with the "right" ionization state, due to the higher dis- 
tance. This result shows in a simple way that, due to the dif- 
ferent dependence of the the physical conditions of the wind 
elements on the initial parameters, a higher luminosity does 
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R(R S ) 

Fig. 5. Ratio between the gas final velocity and its escape ve- 
locity for several values of the Eddington ratio, cedd- The red 
dashed line reproduces the same profile for the "case study" 
shown in Fig. 3. 



not automatically imply a faster wind. In our case the intu- 
itive increase of the wind velocity is observed increasing the 
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Eddington ratio from 0.2 to 0.5, but not with the further in- 
crease to 0.9. 

Density of the gas (Fig. [6j: At small values of the density 
(left plot, n-i = 3, to be compared with nq = 20 in the stan- 
dard case discussed in the previous Section), the gas is ove- 
rionized even at large radii, and it is therefore impossible to 
launch a wind. Increasing the density (nj = 10, right plot) 
the ionization parameter decreases, making it easier to push 
the gas effectively through the line-absorption driven force. 
An even higher density shifts the wind towards inner radii, 
because the over-ionization region is limited to a smaller re- 
gion. 

The ratio fx between X-ray and disc radiation (Fig. [7j: Too 
high a value of fx (in our example, fx =25%, corresponding 
to a 2 keV to 2500A slope of aox = 1-3) causes a higher 
ionization, and make it more difficult to accelerate the gas up 
to the escape velocity. Smaller values (fx =5%, fx =10%, 
corresponding to a x = 1-6 and a x = 1-5, respectively) 
mean a smaller ionizing continuum, and therefore a lower 
ionization factor for the inner gas of the wind, which can be 
accelerated up to the escape velocity. This effect was first 
noted by Murray & Chiang (1995). 

Black hole mass (Fig. [8]): The wind properties strongly de- 
pend on the black hole mass. The dependence is quite com- 
plex, since the different physical parameters (ionization state, 
UV flux, disc temperature) scale in different ways with the 
black hole mass. Varying the black hole mass and leaving the 
other parameters as in our case study, we see that the higher 
the mass, the closer in is the wind. For the highest mass 
value (Mbh = 10 9 M ) no wind is launched with the adopted 
choice of the initial parameters. The main physical driver of 
this behaviour is the ionization parameter, U oc L/R 2 . Since 
L oc M B h (for a fixed e E DD), and the distance in physi- 
cal units is R «M BH (for a fixed value in units of R s ), the 
ionization parameter at a given distance in units of R s de- 
creases with M B h- The distance (in units of R s ) at which 
we have the "right" balance between ionization state and 
UV irradiation therefore decreases with increasing black 
hole mass, reaching values lower than our inner radius 
for M B h = 10 9 M . In this case, it may be possible to ob- 
tain a wind from radii R<R W . This scenario will be stud- 
ied in a wider analysis of the parameter space which will 
be presented in a forthcoming paper. Here we only note 
that this result does not imply that a wind is impossible 
at high masses: an example of wind from a 10 9 M Q black 
hole is presented in the next Subsection, where we show 
that a different choice of the initial density leads to higher 
ionization parameters, shifting the "wind zone" towards 
higher radii. In particular, we will show that an initial 
value of the density Nq=\ (i.e. 20 times smaller than our 
baseline case) implies that a wind can be launched only 
with black hole masses of the order of 10 9 M . 



3.5. Parameter survey 

In order to test the ability of QWIND to perform complete sur- 
veys of the parameter space, we run a grid of 11x10 mod- 
els, adopting the parameters of our baseline study, and varying 
the black hole mass and the Eddington ratio in the range 10 7 - 
10 9 M Q and 0.1-1, respectively. The result is shown in Fig. [9] 
The allowed region for the wind launching extends down to 
£edd ~ 0.2, showing that a radiation wind is possible with lu- 
minosity much lower than the Eddington limit (as also shown in 
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Fig. 6. Ratio between the gas final velocity and its escape 
locity for several values of the initial gas density, rtj^. The 
dashed line reproduces the same profile for the "case stu 
shown in Fig. 3. 



more detail for a single case in Fig. [5J5). Clearly, in these c£ 
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Fig. 7. Ratio between the gas final velocity and its escape ve- 
locity for several values of the ratio fx between the X-ray and 
UV flux. The red dashed line reproduces the same profile for the 
"case study" shown in Fig. 3. 



Fig. 8. Ratio between the gas final velocity and its escape veloc- 
ity for several values of the central black hole mass, Mbh- The 
red dashed line reproduces the same profile for the "case study" 
shown in Fig. 3. 



the enhancement of the radiation force due to line absorption As a further example of the application of QWIND we 
plays a fundamental role. also tested the scenario of a hot, low density wind with a less 
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Fig. 9. Results of a survey of the cedd-Mbh plane, with the other 
initial parameters equal to those of the baseline model. The "+" 
indicate that a wind is launched. 



substantial mass outflow. To test this, we run a QWIND grid 
analogous to the one described above, but with an initial density 
20 times smaller (% = 10 7 cm -3 ). The results are shown 
in Fig. [10] In general, a too low density implies a too high 
ionization, and the wind cannot be launched. The situation is 
different only at very high masses, where the distances from the 
central source are large enough to provide the requested balance 
between ionization parameter and intensity of the UV source. 
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Fig. 10. Results of a survey of the £edd - Mbh plane, with an 
initial gas density rim - 10 7 cirT 3 , 1/20 of the baseline value. 
The "+" indicate that a wind is launched. 



4. Analytic approximations 

Since our main aim is to understand which physical conditions 
allow wind acceleration, we can search for further simplifica- 
tions which allow an even simpler treatment, based on an ana- 
lytic model. 

In this Section we will make three extra approximations to 
the dynamical equations of the wind, which allow the equation 
of motion to be integrated analytically, so that the final velocity 
can be calculated directly. The approximations are: 

1. an energy source with a constant inclination with respect to 
the wind. For the X-rays, this is the same as assumed in 
QWIND (spherically symmetric emission). For the UV, this 
is equivalent to the assumption of a fixed inclination angle of 
the wind stream lines. 

2. neglect the gravitational term in the dynamic equation, as 
once the line-driven acceleration mechanism becomes effec- 
tive, it likely dominates the gravitational term by a factor of 
10 or more (see Section 4.3). 

3. line driving is well reproduced by the simple equation 
M(t) — K x r a , which is valid only when the wind is not 
overionized, but we require a low ionization parameter as a 
condition for the acceleration of the wind (see below). 

Our aim is not to study the details of the wind solutions, 
but rather to explore the parameter space in order to understand 
which are the initial conditions required in order to have acceler- 
ations up to several thousand km s _1 . The analytic treatment lets 
us search for the existence of wind solutions by finding regions 
of parameter space that satisfy two conditions: (1) that the final 
velocity exceeds the escape velocity; and (2) that the gas must 
not be overionized. 

Each of these three approximations makes wind solutions 
more likely: (1) the spherical approximation implies that the UV 
radiation is not decreased by the geometrical projection factor of 
the disc, (2) neglecting the gravitational term obviously favors 
fast outflows, and (3) the force multiplier assumption maximizes 
the effect of line-driving. As a consequence we will obtain wind 
solutions in larger regions of the parameter space than found us- 
ing the more precise numerical approach described in Section 3. 
We will show that in practice the analytically allowed param- 
eter space is close enough to the results of the numerical code 
to make the analytical results useful in delimiting large regions 
of the multi-dimensional initial parameter space where it is not 
possible to accelerate a wind. 

We discuss the third assumption, the conditions on the ion- 
ization state, in the next Subsection. The integration of the equa- 
tion of motion, and condition on the final velocity are discussed 
in Section 4.2. We then show the results for a limited set of ini- 
tial parameters, and finally we compare our results with those 
obtained with QWIND. 

4.1. Relation between M(f) and £ 

Here we show that the equation of the force multiplier, M(t) can 
be simplified into two regimes, depending on the ionization fac- 
tor: £, < 100 : M(f) = Kr a ,£, > 100 : M(t) = 0. We then discuss 
the equation describing the motion of a gas element inside the 
wind. 

In order to estimate the range of values of t for which the 
correction factor of M(t) becomes important, and therefore M(t) 
saturates to its asymptotic limit, we expand Eq. 3 for small f, 
obtaining: 

M(t) ~K(l- a)r] a MAX { 1 - a7] MAX t) (5) 
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This shows that in order to have M(t) 
optical depth must be 



t « 



1 



aVMAX 



the effective 



(6) 



From Fig. [T], for £ < 100 this implies t< 10~ 6 . This condition 
can also be easily obtained from Fig. [2] For an outflowing wind 
reaching a radial velocity of ~ 10 4 km s _1 in a distance of 10 16 - 
10 17 cm, typical AGN values, we have an average ^ ~ 10~ 7 - 
10~ 8 s" 1 and so, from Eq. 1, t ~ 7 X 10~ 3 r£- 5 n 8 , where T = 
10 6 T(, K and n = 10 8 «8 cm -3 . It is easily seen from Fig.Q]that 
if the ionization factor £ is lower than 100, M(t) is never close 
to the saturation regime for reasonable values of the density and 
temperature, and M(t) = K x r a is a good approximation. K is 
constant within a factor ~2 in the range of interest (Fig. 1). 

If £, > 100 the situation is completely different since it is 
obvious from Fig. Q] that the line absorption becomes negligi- 
ble compared with electron scattering and in this condition no 
wind can be launched unless the luminosity is super-Eddington. 
This is why the wind models of Murray & Chiang (1997) re- 
quire additional, "hitchhiking " gas, to shield the gas that was to 
be accelerated, reducing Therefore, in the following we shall 
assume f < 100 in studying the dependence of the wind from 
the other parameters. 

Now let us consider a gas element in the outflow at radial 
distance R from the center. The gas element is subject to grav- 
itational force, Fq = GMBHtn p R~ 2 , and to the radiation force, 
F r = L(AnR 2 cy 1 cr T (M(t) + l). The force multiplier M(f) is given 
by Eq. 3, and depends on the density profile of the gas making 
up the wind. We assume a profile n,-(/?) = n/jy X (R/Rm)^> as m 
our numeric code. 

The flux at a given gas element inside the wind is given by 
FLo/(4nR 2 )xA(R), where Lq is the intrinsic disc luminosity, F is 
the geometrical factor due the disc inclination (which is assumed 
to be constant in this analytical treatment) and A(R) takes into 
account the absorption between the gas element under consid- 
eration and the luminosity source. While a detailed treatment of 
the absorption term would require a complete solution of hydro- 
dynamical and radiative transfer equations, we can make some 
useful approximations: 

- that the inner part of the gas is able to reach high enough 
above the disc in order to shield the gas farther out. This 
assumption is supported by the results of QWIND (§3), and 
has the same effect as the assumption of hitchhiking gas by 
Murray & Chiang (1995). 

- that each ring of gas with thickness AR contributes to radi- 
ation absorption with a column density of ~ nAR. This is 
obviously a lower limit, since the wind is not perpendicular 
to the central radiation. However, the correction is small and 
can be neglected. 

With these two assumptions, the absorption term can be writ- 
ten as: 



A(R) = e 



(7) 



where cr T is the Thomson cross section. Assuming the den- 
sity profile as in our numeric code, n(R) = x (R/Rj^)~^, we 
have: 



We can now write the equation of motion for the gas, simply 
requiring that the force on a gas element is equal to the differ- 
ence between the external radiation force and the gravitational 
force. As in the QWIND code (§3), we neglect the internal pres- 
sure term, since we want to discuss the motion of the ion in the 
supersonic part of the wind, where gas pressure cannot affect the 
dynamics significantly. Furthermore, we assume that only radial 
forces are present. This implies the conservation of angular mo- 
mentum, and a decoupling of the radial and tangential equations. 
The radial equation of motion can then be written, putting €edd 
for the Eddington ratio and R;n for the initial radius, as: 



Vr 



dVR 

dR 



<S>RAn + <t>, 



CRAY 



(9) 



where the left hand side gives the radial acceleration, while 
the right hand side is made up of two terms: 

1. the line-enhanced radiation force, given by the flux multi- 
plied by the force multiplier M(t), and by a geometrical fac- 
tor, F((f>), dependent on the inclination angle <p of the disc 
with respect to the gas element: 



®rad = ^-£EDDA(R)F(cf>)K(nv TH( T T \A\- 1 r a 
R dR 



(10) 



2. the gravitational force, corrected by a factor (1 — ■#), obtained 
from angular momentum conservation, assuming that the gas 
element moves with the Keplerian velocity at the initial ra- 
dius Rj, and decreased by the continuous-radiation force, due 
to Thomson scattering): 



®crav = -^(1 - |)(1 - e EDD F{4,)A(R)) 



(ID 



A further fundamental condition is given by the mass con- 
servation condition, M - AnR 2 nv = const, which implies that 
the only stable solutions are those through the critical point, as 
first discussed by CAK75 and, in a context similar to ours, by 
Murray & Chiang (1995). (See also Lamers & Cassinelli, 1999, 
for a textbook treatment.) 

In general, Eq.QT]can on ly t> e solved using numerical codes. 
However, in our case we are interested in some special situations. 
In particular, the wind solution we are seeking is characterized 
by a sudden and strong radial acceleration (as demonstrated by 
the results of QWIND shown in Section 3), which changes the 
gas motion from slow and quasi-vertical to radial and fast (with 
velocity higher than the escape velocity). We can assume that in 
this accelerating phase the external force is dominant with re- 
spect to the gravitational term, which can be neglected. We also 
require that the final wind velocity, \>f, exceeds the escape ve- 
locity, \'esc = V2GM/? -1 , within a distance from the center of 
a few 10 3 Rs- We will then check whether our solutions ful- 
fill the conditions described above. This consistency check will 
only show whether our wind solutions are acceptable, but other 
solutions are also possible. For example, we expect that for a 
large set of initial conditions, the gas, after initially rising will 
fall back onto the disc (in these cases the gravitational term is 
obviously not negligible). 

Neglecting the gravitational term, and requiring mass con- 
servation, we obtain: 



R 2 v R d -^ 
dR 



A(R) = e 



-§§ „ 8 fl 100 M 7 [l-(^)'- 



(8) 



where 

tj = (vTHCTTV i n i R 2 r a KF(<p)A(R)GMe EDD 



(12) 



(13) 
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Since we are studying solutions with rapid radial accelerations, 
which quickly make the gas motion radial, we can assume that 
the geometrical factor is constant (F((p) = Fo). Furthermore, 
we neglect the dependence of A(R) on the distance R, assum- 
ing A(R)=A(1.2R,). This approximation is reasonable, con- 
sidering that our simulations with QWIND (Section 3) show 
that when a wind is successfully launched, the escape veloc- 
ity is reached at a distance from the center R ES c only slightly 
larger (by on average ~20%, from Figs.5-8) than the initial 
distance of the stream line, /?,-. With these assumptions, Eq.[T2l 
can be easily integrated between an initial radius, /?,-, and a final 
radius, Rf. The condition v E > v E can then be written as: 

n^n'n-^-] >gmr f 1 (14) 

R F 

Finally, we adopt the initial density profile discussed above, and 
we assume Rf-2Rj. A choice of Rf is needed in our treatment, 
because having neglected the gravitational term implies that the 
above equation is always satisfied for Rf — > oo. We assume 
Rf-2Ri, which, based on the above considerations on R E sc, is 
a rather conservative assumption (in all the cases analyzed with 
QWIND we found R E sc<2Ri), in the sense of allowing a larger 
parameter space for wind launching than in our analytic treat- 
ment. This is in agreement with our general approach, as dis- 
cussed at the beginning of the present Section. 
We can then rewrite the above equation as: 

[v TH cr T v I n I R 2 i {GMy l )Y a kF AeEDD > 1 (15) 

Finally, adopting a proper parametrization for the physical 
quantities and substituting the numerical values for the con- 
stants, and assuming a=0.6, we obtain: 

e EDD > 1.1 x \Q^A- l F^[^JT lVim Rl m) M- l f 6 (16) 

This condition, when combined with the requirement that the 
gas not be over-ionized (§4.2), will allow us to put constraints 
for the existence of wind-like solutions in a black hole mass - 
Eddington ratio plane, as we will show in Section 4.3. 

4.2. Effects of X-ray absorption 

If the inner part of the wind is too close to the X-ray source, an 
over-ionization problem arises. An approximate treatment can 
be obtained by considering the X-ray cross-section of the wind 
gas as in POO. Using the approximation described in Section 
3, a fundamental requirement for a radiation-driven wind to be 
launched is that £=L/(n(R)R 2 )<10 5 , where L is the ionizing 
radiation. Adopting the same parametrizations as in the pre- 
vious Subsection, the above condition can be written as: 

1.5 x 10 & f x n, l R-^M- l A(R)e EDD < 10 s (17) 

where A{R) is the fraction of non-absorbed X-ray radiation, and 
is given by Eq[8] the cross section being the same (for £ > 10 s ) 
as for the UV disc radiation. f x is the fraction of the bolomet- 
ric luminosity emitted in the X-rays. In the work of Stevens & 
Kallman (1986) the objects considered as sources of radiation- 
driven winds were X-ray binaries, with a rather hard spectrum. 
As a consequence, only photons with E>1 keV were taken into 
account in their estimate of the ionization parameter. 

In our case, the typical intrinsic X-ray spectrum of an AGN 
is dominated by soft emission in the 0.1-1 keV band, which gives 
a significant -if not dominant- contribution to the ionization pa- 
rameter. Therefore, we calculate the ionization parameter using 



the luminosity above 0. 1 keV. As a standard value we assume 
that of Elvis et al. (1994), where the 0.1-100 keV luminosity is 
~ 15% of the total (i.e. optical-UV-X-ray) emission of an AGN. 

Assuming the density profile of the SS disc n(R) = x 
(R/Rin)~P equation [17] can be finally rewritten as: 

e EDD < 0.0044(-^-)- I M 7 « 8 (^-)^A(/?)- 1 ^ 00 (18) 



4.3. Analytical Results 

The constraints from EqQ~6]and Eq [18] delimit the region of the 
black hole mass, Eddington ratio (M BH , € E dd) parameter space 
in which an outflowing wind can arise. 

We can use EqsQ~6] and [18] to plot the allowed wind- 
launching regions in the cedd - Mbh plane. Since Eq. [18] de- 
pends on the initial radius of each streamline, for a given set of 
initial parameters (wind inner radius R/n, density at the inner ra- 
dius R[n, gas temperature T and X-ray to optical/UV ratio fx) 
we compute the allowed region for each single streamline, and 
then we plot the convolution of all these regions. In this way, 
the meaning of the "analytic parameter surveys" are analogous 
to those plotted in Figs. |9l and [TOl 

We assume a baseline set of initial parameters (Table 13. 2K 
as used in EqsOS] and ED fx =15%, T = 2 x 10 6 K, n IN = 
10 8 cirT 3 at Rin - 100 Rs- These parameter values are those 
required for the gas responsible for the observed UV and X-ray 
'Warm absorber' absorption lines (Nicastro et al. 1999, Netzer 
et al. 2002, Rrongold et al. 2003) in quasar spectra. They are 
also the physical conditions needed for a gas confining the BEL 
clouds (which have n ~ 10 10 cm -3 and T ~ 10 4 K, Osterbrock 
1989) as noted by Turner et al. (1994) and Elvis (2000). 
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Fig. 11. Results of our numerical analysis for our baseline set of 
parameters. The green (shaded) region indicates where a wind 
can be launched. The analytic estimate is compared with the nu- 
merical result obtained with a QWIND survey. 

Figure QT] shows the allowed {Mbh, £edd) parameter space 
for the baseline parameters, superimposed to the results of 
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QWIND for the same set of parameters (already shown in Fig. [9] 
The two interesting results emerging from this plot are: 

1) The allowed region based on our analytic approximation is 
larger than that obtained with the more detailed approach based 
on QWIND. This confirms that our approximations do not miss 
any possible wind solution, in agreement with our expectations 
(Section 4.1). 

2) The results of the analytic approximations are useful to ex- 
clude a significant part of the parameter space where a wind 
cannot exist. This approach can then be used for a genaral explo- 
ration of the parameter space, leaving the more time-consuming 
numerical approach for a smaller set of parameters. This was not 
obvious a priori, since a too strong approximation could have 
lead to a too large (and therefore useless) allowed region. 

While a complete investigation of the parameter space is not 
the purpose of the present work, it is useful to show the results 
of the analytic approach for a small set of initial parameters, as 
examples of the possible uses of this technique. 

We show in Fig. [T2]the allowed regions in four cases, ob- 
tained from the baseline case changing one or two initial param- 
eters, in order to test four different physical situations: 
(a) BELR wind: It is interesting to investigate whether BEL gas 
could be accelerated directly by the radiation force. Fig. [T2"k 
shows that BEL gas cannot be accelerated to high velocities by 
the radiation force, except for relatively low mass black holes 
(Mi < 10) and high Eddington ratios (excpet for the lowest 
masses). The physical reason is that in an overly dense gas 
the Doppler shift due to the radial acceleration is insufficient 
to prevent the gas from being self-shielded against line absorp- 
tion. Thus in a too dense gas the effective optical depth, t, is 
never small enough to make the force multiplier M(f) much 
higher than 1, for reasonable values of the radial velocity gra- 
dient (^ > 1CT 6 s _1 ). We stress that this result does not imply 
that a wind containing a cold and dense phase is possible only at 
low masses. Indeed, if the BEL gas is confined by a hot wind, it 
is expected that the pressure of the warm gas would drive them 
to similar velocities. However, a more detailed study is needed 
to test this statement. Here we only note that in a wind like the 
one described in Elvis (2000) most of the kinetic energy is in 
the warm phase (80-90%). Therefore, from an energetic point 
of view it is likely that the dynamics of the whole wind (cold 
phase, i.e. the BEL gas, and warm phase) is determined by the 
warm phase. 

(b) Hot wind: Fig. [T2b shows the allowed region for a hot, 
low density wind. Only at the highest masses can such a wind 
be launched, because of the over-ionization due the low wind 
density. The ionization parameter of the wind decreases at high 
BH masses because in our scheme the distance of the stream 
lines is always the same in units of R$ . Therefore, increasing the 
BH mass implies a linear increase of the luminosity (at a given 
£edd) and a quadratic increase of the physical distance, resulting 
in a linear decrease of the ionization parameter with Mbh- 

We note that this case is the same as in the second QWIND 
survey (Fig. 10). Again, the allowed region obtained with the 
analytic approach is larger than that predicted with QWIND 
(Fig. HE 

(c) X-ray weak source: The fraction of the bolometric lu- 
minosity emitted in the X-rays, fx, is important in the deter- 
mination of the allowed region for a wind. The value used in 
the standard parameters set, fx-15%, is in agreement with the 
observed Spectral Energy Distributions (SED) of PG quasars 
(Elvis et al. 1994, Laor et al. 1997). However fx depends on 
luminosity (Zamorani et al. 1981, Yuan et al. 1998, Steffen 
et al. 2006) and so higher values, ~ 20 - 25% are found in 



nearby Seyfert galaxies. High values are also found in radio loud 
quasars (in this case however beaming could play a role). On 
the other side, lower values, ~5% are found in high luminos- 
ity (L 250oA > 10 31 erg s _1 Hz -1 , Yuan et al. 1998) and high 
redshift (z>4, Vignali et al. 2002, Steffen et al.2006) quasars. 
In Fig. [T2b we show the case for an X-ray quiet source, with 
fx = 5% and the other parameters as in the baseline study. The 
allowed region is increased, because lower X-ray emission im- 
plies a lower ionization parameter. This is important especially 
at low BH masses, where the requirements to avoid gas over- 
ionization are more stringent. 

(d) Wind Inner Radius: We changed the inner wind radius 
from 100/?s to 3QQRs (Fig. [TZtl). The allowed wind region is 
smaller than in our baseline study, because part of the parameter 
space is excluded both at high BH masses (as is apparent from a 
comparison with Fig.fTTI). where a wind cannot be launched be- 
cause the physical distance of the gas is too high for the radiation 
pressure to be effective. 

Summarizing, our results indicate that only a limited range 
of physical conditions can give rise to a radiation driven wind in 
AGNs, but also that these conditions do exist, and are the same 
as those observed in AGN warm absorbers, and needed for a 
warm gas that pressure confines the BEL gas. 

Furthermore, we showed how our simple approach can pro- 
vide an immediate physical understanding of the dependence of 
the wind on the individual initial parameters. 

The analytical results presented here have several important 
limitations, as noted at the start of Section 4. In addition, our 
treatment is based on a condition of existence for the solutions 
of the equation of motion of a wind element. This is a less strin- 
gent condition than requiring an actual wind solution for a given 
set of initial parameters. As a result, real wind solutions are pos- 
sible for smaller parameter regions than those shown in Figs. QT| 
and [12] Importantly though, we do not expect any real solution 
outside these regions. The exception is for £edd ^ L for which 
a wind can be accelerated through electron scattering alone. 

5. Comparison with previous results 

A comparison of our results with those obtained with more com- 
plex hydrodynamical codes is useful to test the reliability of 
QWIND. In particular, we refer to the work of POO. in which 
a radiation-driven wind from a SS disc is studied, solving the 
wind equations through a 2-dimensional hydrodynamical code. 

The properties of the solutions are remarkably similar. In 
particular: 

- The velocity profile in POO (their Fig. 2) is remarkably sim- 
ilar after the wind has bent towards the radial direction. In 
particular, when the external force becomes effective, and 
the gas changes its direction, the velocity very rapidly in- 
creases up to its maximum value (around 10 4 km s~'). This 
is in agreement with our main approximation, i.e. the sepa- 
ration between a internal pressure dominated phase and an 
external radiation force dominated phase. 

- The ionization parameter profiles are quite similar. This is 
not surprising, and indeed is only a check of self-consistency, 
since we have adopted the same assumptions for the X-ray 
cross-section as POO. As discussed above, in this scenario 
the ionization parameters rapidly drops from large (> 10 5 ) 
to very small (< 10 3 ) values across a thin layer of Nh ~ 
10 22 cm" 2 . 

- The final values for the ejected mass are comparable. We find 
for our "case study" a value of ~ 1 M yr -1 close to that of 
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Fig. 12. Results of our numerical analysis for several choices of the initial parameters. The red lines show the total luminosity of the 
source. The green (shaded) regions are those allowed for wind launching. Upper left panel: BELR case; upper right panel: hot wind 
case; lower left panel: wind from an X-ray quiet source; lower right panel: wind with a large inner radius (R^ = 300 Rc)- 



POO for the same black hole mass (we note that the mass loss 
rate scales as M\ H ). 



Summarizing, the comparison with the results of POO shows 
that our model is in qualitative agreement with a much more 
complex and more precise treatment. Given the huge uncertain- 
ties in the initial conditions (both regarding the structure of the 
disc and the mechanism of initial launching of the gas) we be- 
lieve our approach is well supported. We remind the reader how- 
ever that our simplified approach should be used as a way to 
quickly explore and physically understand the possible solu- 
tions, while a more detailed treatment is needed if one is inter- 
ested in the precise values of the parameters involved or in more 
detailed structure. 



6. Conclusions, and Future Work 

We have shown that a radiation-driven wind can be accelerated 
to velocities of about 10 4 km s _1 from a Shakura-Sunyaev ac- 
cretion disc at a distance from the central black hole of the order 
of that of the BEL clouds, for densities and temperatures in the 
range observed in quasars warm absorbers. 

Thanks to our simple approach, we were able to make an an- 
alytic estimate of the parameter space in which a wind solution is 
allowed, and to build a fast code, QWIND, which can be used to 
better investigate the dependence on the several unknown initial 
parameters. 

The main limitation of our study is the assumption of a sub- 
stantial initial velocity (a few 100 km s _1 ) which makes the wind 
supersonic from the beginning. We did not discuss the physical 
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mechanism through which the disc provides the initial kinetic 
energy to the gas. 

Our main results are: 

1. Our model reproduces the global results of more complex 
hydrodynamical codes. This makes us confident that our 
approach can be used to explore the huge space of initial 
conditions, in order to understand in which cases a wind 
can be launched. The results of our analytical study are in 
agreement with those of the numerical code QWIND, in the 
sense that parameters space regions excluded by the analyti- 
cal analysis are also excluded with the more precise numeri- 
cal approach. 

2. For initial parameters typical of AGN warm absorbers, a 
wind can arise only for relatively high accretion efficiencies 
(cedd > 0.3) and in a narrow range of disc radii. This is 
due to overionization at too small radii, and to a too small 
UV radiation and initial velocity at too large radii, and is an 
inherent feature of radiation driven disc winds in AGNs. 

3. The inner part of the gas, which is too overionized to be ef- 
fectively accelerated, forms a "failed wind" which removes 
the X-ray radiation from directions close to the disc plane, 
allowing a fast decrease of the ionization parameter with 
the distance form the central source, and so allowing the 
line driven wind acceleration at larger radii. This compo- 
nent is analogous to the "hitchhiking gas" of Murray & 
Chiang (1995). 

4. We have shown how QWIND can be used to gain a simple 
physical understanding of the dependence of the wind prop- 
erties on the initial parameters, such as the X-ray to UV ra- 
tio, the Eddington ratio, the initial density and temperature, 
the central black hole mass. Some of the results are at first 
counter-intuitive, as, for example, the non-monotonic depen- 
dence of the wind final velocity on the Eddington ratio: in 
some cases more luminous sources can have slower winds, 
with the other physical parameters being held constant. 

5. The terminal velocity of the wind is typically of the order 
of 1 - 2 10 4 km s _1 , and the cone covering factor is a few 
percent of the lines of sight. This is in reasonable agreement 
with observations of BAL quasars. 

6. The wind angle above the disc is substantial (~ 20deg), and 
so is neither equatorial nor polar. Some ~ 35% of the view- 
ing directions will see the X-ray and UV continuum source 
through the wind. 

Our findings could have deep implications for wind- 
embedded BEL models, like that of Elvis (2000). Black holes 
accreting in conditions outside the narrow allowed ranges would 
not have winds, and so neither absorption lines, nor the BELs 
would be observed. Such 'naked quasars' would be picked up 
from their non-stellar colors in the Sloan Digital Sky Survey 
or by their high X-ray flux in X-ray surveys. A few "lineless" 
quasars have been found but they are rare (e.g. McDowell et 
al. 1995). If AGN winds are radiation-driven, then some self- 
regulation mechanism must be involved in the accretion process, 
to keep the AGNs at the conditions needed for a radiation driven 
wind. 

Our main aims at this stage were to check the reliability of 
our approach, and to give significant examples of its use. The 
next major step, which will be presented in a forthcoming paper, 
is a systematic survey of initial parameters, to give a more gen- 
eral view of the physical conditions needed to accelerate winds 
by line driving. In particular, we will investigate how the results 
are affected by changes in the wind inner radius and in the initial 
density and velocity profiles. 
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